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Abstract. 

A parametrically forced Bose-Einstein condensate (BEC) is studied in the mean 
field approach for the case of a general moving optical lattice. The interaction between 
the atoms in the condensate and the time dependent lattice potential leads to a 
novel propagating superfluid matter wave, which can be controlled through chirp 
management. This system, when placed in a trap, accelerates and undergoes rapid 
nonlinear compression, controlled by the chirp. The density achieves its maximum, 
precisely when the matter wave changes direction. A dynamical phase transition is 
identified, which takes the superfluid phase to an insulating state. The exact expression 
for energy is obtained and analyzed in detail to gain physical understanding of the chirp 
management of the sinusoidal excitations and also the dynamical phase transition. 
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1. Introduction 

Bose-Einstein condensate (BEC) in a periodic potential is an area of active research 
[USEIIUEIEIITIIHIIHI ^ ■ ^ * s an interdisciplinary field, which has connections with 
many areas of physics: electrons in crystal lattices, polarons and photons in optical 
fibers, gauge theories and exotic phase transitions [H] [121 EB], to mention a few. The 
nature of the phase transition is different for deep and shallow optical lattices (OL). 
In the former metal to insulator quantum phase transition [12j, driven by the 

quantum fluctuations can occur, whereas, in the latter case, the same can occur through 
a classical dynamical phase transition [HJ [151 EE EE] • The advantage of studying the 
condensate in optical lattice is that one can control and tune all the relevant parameters, 
such as lattice depth, its geometry and the interaction strength between the atoms 
PS [19j EE] EH1 121] ■ Optical lattices (OL) can be made time dependent by temporal 
variations of the depth and the width of the lattice [22J. The lattice can be moved at 
a velocity v = (X/2)5u, where A/2 is the spatial period of the periodic potential and 
8v is the stable detuning between the two laser beams, which also allows the lattice to 
accelerate. The amplitude of the lattice potential can be tuned through the following 
relation: V$ ~ 2HHr/ Er, £Ir being the Rabi frequency and Er = 2h 2 7i 2 jm\ 2 is the 
recoil energy. 

A number of recent studies have investigated the condensate dynamics, in presence 
of time modulated optical lattices [23] [2U [25] [26] . Madison et al , studied quantum 
transport of atoms in presence of an accelerating lattice [23]. They observed resonances, 
when a small oscillatory trap is added to the system. The number of atoms decreases 
rapidly, when the resonances occur. The spatiotemporal dynamics of the condensate 
has been investigated in presence of moving optical lattice, using direct perturbation 
and Melnikov-function method [25]. The existence of sub diffract ive solitons is shown in 
[26] , for a time dependent shallow optical lattice. These authors also showed that the 
variation of the chemical potential can lead to a chaotic system, whereas, by controlling 
the intensity of the lattice potential, the chaos can be suppressed. It is worth pointing 
out that, Fallani et al , have observed dynamical instability in a one dimensional moving 
optical lattice [22]. 

We present here a detailed study of the dynamics of BEC in presence of a general 
time dependent optical lattice with both attractive and repulsive atom-atom interaction. 
The exact sinusoidal solutions are naturally chirped, which allows them to be compressed 
and accelerated in a controlled manner. We further consider the BEC in the combined 
presence of optical lattice and a harmonic oscillator (HO) trap. The presence of the 
trap does not affect the form of the solutions, but drastically changes the dynamics of 
the condensate. It is found that the BEC can accelerate and undergo rapid nonlinear 
compression, when the matter wave changes its direction. The density of the condensate 
is found to be maximum at this point. In order to gain a deeper insight in to the 
response of the condensate to various time varying parameters, we explicitly compute 
the energy of the system. In the combined presence of a harmonic trap (HO) and optical 



Chirp Control of Sinusoidal Lattice Modes in Bose-Einstein Condensate 



3 



lattice, BEC shows a resonant increase in energy at certain points in the scaled time 
variable. This is akin to a recently observed resonant behavior in an optical lattice, 
where energy transfer takes place between two bands [27] . In case of the expulsive 
oscillator, the energy initially increases with time and takes a maximum value, after 
which it decreases rapidly. In case of the moving lattice, as the BEC expands, the 
energy monotonically decreases with time. It is also shown that this system exhibit 
dynamical superfluid insulator transition (DSIT), where superfluidity breaks down and 
the condensate transits to an insulating phase. The stability of the obtained solutions 
are carried out using Vakitov-Kolokolov criterion (VKC) [281 [23 EQ], where we found 
that solutions are stable for attractive atom-atom interaction. 



2. Sinusoidal excitation in time dependent optical lattice 

In presence of a shallow optical lattice, the mean field Gross-Pitaevskii (GP) equation 
well captures the condensate dynamics: 

l% ^t = ("2^ V2 + Vext + U ^~ p (*)) *■ C 2 - 1 ) 

Here, V ext = Vno(x,y) + V\(z,t) is the external trapping potential, U is the strength of 
the atom-atom interaction and u(t) is the time dependent chemical potential. To have 
a cigar shaped BEC, one applies a strong oscillator trapping potential of frequency 
UJ ±, Vho( x iD) — \rnujx\x 1 + y 2 ), along the transverse direction. Assuming tight 
transverse confinement, the trial wavefunction is taken as, ^ = ip(z,t)(f) Q (x,y), where 
4>o is represented by a Gaussian ansatz [3f l 132 1 133]: 



x 2 +y 2 



(2.2) 



with a± = sj^. In the weak coupling regime, one finds that the quasi one dimensional 
GP equation takes the form [32 [351 HE S3 EE] : 



^ = + V x {z, W + g\^ - u(t)1>. (2.3) 

Here, V\(z,t) = Vol(z, t) = Vo(t) cos 2 (T) is the time dependent optical lattice potential, 
with T = A(t)z and u(t) = fiA 2 (t). The width A(t) can be controlled by the frequency 
of the laser beam. Vo(t) is the scaled time dependent amplitude of the lattice potential, 
which can be tuned by the intensity of the laser beam. The reduced two-body interaction 
is given by, g = ^^(^.)U, where, k = 2n/\. The time, spatial coordinate and the 
wavefunction are, respectively, scaled as t — > m/(hk 2 )t, z — > z/k and ip — > y/ktjj. The 
chemical potential v(t) and the amplitude of the lattice potential have been normalized 
in terms of the recoil energy Er = ^p-. 

We assume the following ansatz for the BEC profile: 

V>0M) = V ' A{t)a{Ty^ z V +4 '^ . (2.4) 

where, x{ z i t) is a nontrivial phase, controlling the supercurrent of the condensate and 
the chirped phase (f>(z,t), is of the form: 4>(z,t) = —\c{t)z 2 . 
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Current conservation yields, |^ = where 5 = ^y(l — 2/i + 2a) (1 — 2/x). The 
chirp parameter c(i) satisfies = c 2 (t); from which all other parameters can be 
determined. The width A(t) is related to c(t) by, 

^ = A(t)c(t). (2.5) 

For the solution to exist, the width and the depth of the lattice are related 
by, Vo(t) = aA 2 (t) and the two-body interaction strength needs to be of the form: 
g(t) = nA(t). The real part of the GP equation, in terms of the density, then reduces 
to: 

1 d 2 a 1 ( da\ 2 , , 9/mN 9 , 

~ 7^7^ -ob - ^ + Ka + acos 2 (T)a 2 + 25 2 = 0. (2.6) 



4 9T 2 

For constant a, one finds a self-similar solution of the form: 

a(T) = a + 5cos 2 (T). (2.7) 

where, a = and b = —-. Notice that, for cr(T) to be positive semi-definite, 

considering the atom-atom interaction to be repulsive, /i > a + 1/2, when a is positive. 
For, negative a, // > 1/2. When the interaction is taken to be attractive, for positive 
and negative a, n > 1/2 and /i > a + 1/2, respectively. 

The periodicity of the density modulation is same as that of the lattice potential 
as seen in figured] Keeping in mind the role of chirp c(t) = l ^° cot , for effective pulse 
compression in optical fiber [39], we consider a time dependent chirp parameter, of 
the same form. Here, A(t) = 1 ^° ot , with the initial condition c(t = 0) = Co and 
A(t = 0) = Aq. The compression of the matter wave is shown in figure la. Chirping 
leads to a rapid compression of the BEC profile and the density tends to a singularity 
when t — > 1/cq. The number of atoms per lattice site, in this case, is given by, 

iV= X(2 ";;" 1) -£s W , (2.8) 

which can be controlled through chirp management. 

Till now, we have restricted ourselves to the static optical lattice, whose width and 
amplitude have been allowed to vary with time. We now allow the lattice to move with 
a velocity v = |L with T = A(t)(x — l(t)). In this scenario, the solution retains its 
original form. The position of the matter wave can be found from the relation: 

/(t)=/ e-^ c (^'=/ (l + Co t), (2.9) 

showing a constant velocity v = IqCq. In the case, with c(t = 0) = —1/cq, we find 
c(i) = — l/(t + Co) and A(t) = A c /(t + Co). As time increases, the width, as well as 
amplitude of the matter wave decreases, leading to the spread of the BEC, as depicted 
in figure lb. 

3. Dynamics in presence of an optical lattice and harmonic trap 



Keeping in mind the fact that most of the atomic systems are inhomogeneous due to 
the presence of the magnetic or optical traps, we now consider the same in a harmonic 




(c) (d) 

Figure 1. Figure la depicts the sinusoidal BEC profile in presence of a static 
optical lattice, where the chirped BEC undergoes rapid compression. Figure lb 
corresponds to a different parameter domain, where BEC undergoes expansion. Figure 
lc represents the density in presence of a harmonic trap, where at the point of nonlinear 
compression, the matter wave changes its direction. Figure Id corresponds to the 
expulsive oscillator, where the BEC spreads. In all these cases, the parameter values 
arc: cq = 1, Aq = 0.5, a = 0.15 and n = 0.8. 



trap, of the form, Vho{ z ^) — ^>M(t)z 2 , in addition to the lattice potential. Control of 
BEC and its solitonic and sinusoidal excitations is an area of significant current interest 
[iQj HD H21 HHJ HU H3 HO]- In a harmonic trap, this is achieved through the time 
dependencies of scattering length, transverse and longitudinal trap frequencies. It has 
been observed that, the center of mass of solitons, as well as its width and amplitude gets 
coupled to the trap parameters, through which they can be accelerated or compressed 
[171 HE], which is useful for coherent atom optics (331 EQl EI]- In the present case, all 
the parameters are controlled by the harmonic oscillator, although the solutions retain 
their self-similar form. The chirp parameter c(t) can be found from the following Ricatti 
equation: 

ct(t) = c 2 (t) + M(t). (3.1) 

We have analyzed the following cases, considering the fact that the oscillator can 
be both regular or expulsive. 

(i) In the first case, we assume M(t) = q 2 , representing a regular trap. Here, two 
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solutions are possible: (i) c(t) = gtangt and (ii) c(t) = — q cot qt. Consequently, the 
other parameters can be obtained. For (i), A(t) = A / cos qt and V (t) = V / cos 2 qt 
and for (ii), A(t) = A / sin qt and Vo(t) = Vo/sin 2 gt, with Vo = a.A\. The 
position of the matter wave is oscillatory: (i) l(t) = locosqt and (ii) l(t) = Iq sin qt, 
respectively, with the initial position Iq. The density profile for this case is depicted 
in figure lc. The sinusoidal wave propagates with a periodicity n and at the point 
of nonlinear compression, the density becomes maximum, wherein, the matter wave 
reverses its direction, as seen in figure lc. Here, the supercurrent takes the form, 

J(t) = mz,t)\ 2 ^ = 2S. (3.2) 

(ii) If one considers M(t) = —q 2 , which corresponds to an expulsive oscillator, the 
Ricatti equation takes the form: c t (t) = c 2 (t) — q 2 . Like the previous case, here 
also we have two solutions: (i) c(t) = — gtanhgt and (ii) c(t) = —qcothqt. The 
width and lattice amplitude for (i) are: A(t) = A / coshqt, Vo(t) = Vo/ cosh 2 qt and 
for (ii) A(t) = A / sinhgt, V (t) = V / sinh 2 qt. The density profile for this case 
is shown in figure Id. As time increases, the amplitude of the propagating wave 
decreases, making the BEC spread out. 

4. Energy of the excitation 

To explore the dynamics of this parametrically driven system, we study the temporal 
behavior of energy, both in the absence and presence of the trap. In figure 2a, the same 
is shown for a static optical lattice, where the energy increases with time for this driven 
system. Figure 2b depicts the behavior of energy for the optical lattice moving with a 
constant velocity v, where the energy is found to decrease rapidly as the BEC expands. 
Interestingly, in presence of a harmonic trap, the BEC shows a resonant increase in 
energy, at some suitable values of the scaled time variable, where it undergoes rapid 
nonlinear compression, as seen in figure 2c. The contribution of the chirping to the 
kinetic energy is solely responsible for this phenomenon. The matter wave changes its 
direction at the point of nonlinear compression. In case of an expulsive oscillator, as 
depicted in figure 3d, one can see that initially energy increases with time and reaches 
its maximum value after which it decreases rapidly. The number density in this case 
turns out to be: 

N = A( 2 -^-^-)-^(sin(2Al) 
V 2k 2k J Atvk K ' 

+ sin(2A(7r -/))). (4.1) 

In case of moving lattice, the number of atoms decreases as a function of time, which is 
similar to the experimental observation in [22] . 

5. Dynamical superfluid-insulator transition 

We now examine the possibility of dynamical phase transition in this system. In 
case of shallow lattice, a dynamical superfluid-insulator transition (DSIT), driven 
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(c) (d) 

Figure 2. Figure 2a shows the rapid increase in energy with time, for a static lattice. 
In figure 2b, the energy decreases with time due to expansion of BEC. A resonant 
increase in energy is observed in figure 2c, where the matter changes its direction. In 
presence of expulsive oscillator, depicted in hgure 2d, energy increases and takes a 
maximum value, after which it decreases rapidly as a function of time. 



by modulational instability, was predicted by Smerzi et al , using the mean field 
discrete nonlinear Schrodinger equation [TJ]. In the following year, Cataliotti et al , 
experimentally observed this classical phase transition in presence of a stationary optical 
lattice [15] . It has been found [IT] that, DSIT can also occur when both two and three- 
body interactions are present. Recently, Fallani et al , have observed this dynamical 
instability in a one dimensional moving optical lattice [22]. Since the dynamical phase 
transition occurs at the point, where the energy of the system becomes non- analytic 
(NA), we have investigated the same by computing the the exact energy for all the 
three cases, separately. The expression for energy is rather too lengthy to be reported 
here, we will only concentrate on the non-analytic term, which for all the three cases is 
given by, 

AA 2 (A 2 a(a + b) - 5 2 ) _-, y/a tan Al . . 

E NA = tan — , (5.1) 

Ja(a + b) Va + b 

In the presence of a harmonic trap, the addition contribution to the energy due to 
the trap arises, is analytic. It is non-analytic at a + b = and a = 0, which respectively, 
lead to 2/i + 2a — 1 = and 2/x — 1 = 0. Both of these give 5 = 0. For the density to 
be a finite positive quantity, we exclude the second point. Taking into account the first 
condition, one can see that the supercurrent vanishes at this point and the superfluid 
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phase transits to an insulating phase: 



foiz, t) = J sin {A{t)z)e l<t> . (5.2) 

V K 

We now check the stability of the above solution by using VKC, which states that, if 
' v > 0, the solutions are unstable and are marginally stable if ^ = 0. Stable solutions 

dN 



correspond to ^ < 0. For all the three possible cases, we found: 
ON A(t) 



(5.3) 



dfj, k 

Therefore, one can note that the stability of the sinusoidal solution depends on the 
chirp, as well as on the nature of the atom-atom interaction. The solutions are stable 
for attractive atom-atom interaction. 



6. Conclusion 



In summary, we presented a detailed study of the parametrically forced GP equation 
in general time dependent optical lattices, with both attractive and repulsive two- 
body interaction. The interaction between atoms and the lattice leads to sinusoidal 
solutions, which can be controlled by chirp management. Both for stationary and the 
time dependent optical lattice scenario, a dynamical phase transition occurs, which 
takes the superfluid phase to the insulating phase. The BEC can be made to accelerate 
and undergo rapid compression by incorporating harmonic trap, apart from the lattice. 
A resonant increase in energy is observed at a certain point of time, when the system 
couples to the harmonic trap. The combined effect of lattice and trap leads to a nonlinear 
compression of the condensate, at certain intervals of time, where the density becomes 
maximum. In case of expulsive trap, the energy increases with time initially and achieves 
its maximum value, after which, it decreases rapidly. We have carried out the stability 
analysis using VKC, and found that the stability of these solutions depends on the chirp, 
as well as on the nature of the interaction. We hope that these observations will find 
experimental realization in the present laboratory setup. 
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